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^..^1 Abstract 

' We calculate a general spectral correlation function of products and ratios of characteristic 

polynomials for a N x N random matrix taken from the chiral Gaussian Unitary Ensemble 
(chGUE). Our derivation is based upon finding a Itzykson-Zuber type integral for matrices 
from the non-compact manifold GI(n,C)/U(l) x ... x U(l) (matrix Macdonald function). The 
correlation function is shown to be always represented in a determinant form generalising the 
known expressions for only positive moments. Finally, we present the asymptotic formula for 
^ the correlation function in the large matrix size limit. 

in 

^ ; 1 Introduction 
lO ■ 

■ A class of random matrices that has attracted a considerable attention recently [Q, |^, ||, |^, ^, |[ 

! 0! H' B ll^l 0' 111' HI' llll ^^"^ so-called chiral GUE, also known as the Laguerre ensemble. The 

corresponding matrices are of the form Z) = ^ ) ' ^'^^'"^ stands for a complex matrix, 

- , with jt being its Hcrmitian conjugate. 

The off-diagonal block structure is characteristic for systems with chiral symmetry. The chi- 
, ral ensemble was introduced to provide a background for calculatin g th e universal part of the 

microscopic level density for the Euclidian QCD Dirac operator, see |15| and references therein. 
Independently and simultaneously it was realised that the same chiral ensemble is describing a new 
' group structure associated with scattering in disordered mesoscopic wires . 

^ I One of the main objects of interest in QCD is the so-called Euclidean partition function used 

to describe a system of quarks characterized by n/ flavors and quark masses my interacting with 
the Yang-Mills gauge fields. At the level of Random Matrix Theory the true partition function is 
replaced by the matrix integral: 

Z„,(M/) - / Pi n det{*i> + m('^)l2^}e-^Try(jt J) 
fe=i 

where V{z) is a suitable potential. Here the integration over complex J serves to mimic the func- 
tional integral over gauge field configurations jist . Then the calculation of the partition function 
amounts to performing the ensemble average of the product of characteristic polynomials of iD over 
the probability density P{J) oc e . To model the simplest case of the sector with zero 

topological charge the matrices J have to be chosen general N x N complex, and the probability 
distribution can be chosen Gaussian as defined by the formula: 

N 

P{J)dJdj'' = const exp (^-NTi (i^^)) H {dReJ^J){dImJ,J) (2) 
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In a more general case of non-zero topological charge the matrices J have to be chosen rectangular 
Jist , see also[|l6|. One may also wish to consider a more general type of potentials We 
do not consider these modifications in the present paper. 

The characteristic feature of the chiral ensemble is the presence of a particular point A = in 
the spectrum, also called the " hard edge" |^ . The eigenvalues of chiral matrices appear in pairs 
±Afe , k = 1, N . Far from the hard edge the statistics of eigenvalues is practically the same as 
for usual GUE matrices without chiral structure, but in the vicinity of the edge eigenvalues behave 
very differently. 

Define the following characteristic polynomials for chiral matrices: 

Z( J, m) = det f ml2N + id) = det ( "'■';f ^ . (3) 

V / \ iJ' mlN / 

To calculate the eigenvalue density p(A) one can use the trace of the resolvent for the matrix 
—iD which is given in terms of the eigenvalues by 

1 ^ 1 
R±{m) = Tr — = 2m^ 



ml2N ± iD fr[ "^^ + K 

The eigenvalue density is then extracted in the usual way as 

p(A) = hm ^Re Tr \ \^^^_^^^ 

In turn, the trace of the resolvent can be most easily found via the identity: Tr(mbl2jv + iD)^^ = 
dmj \ mf=mt ■ Wc therefore see that the interest in spectral properties of chiral matrices 

suggests to consider a general correlation function containing both products and ratios of the 
characteristic polynomials: 

1=1 ^ ' J 

where M/ = diag {rn^f' , ■ ■ • , "^/"^^^ : — diag {^^^ , ■ ■ • , '^i^"''''^ and the average is taken over the 
ensemble Eq.(H). In fact, correlation functions of such a type contain the most detailed information 
about the spectra of random matrices and can be used in various contexts and applications |]l7| . 

There are several analytical techniques for dealing with the correlation function ^Jlfc^i Z{J , m,k) 
of only products of characteristic polynomials. The majority of the results obtained so far resorted 
to exploitation of the orthogonal system of Laguerre polynomials, see [Q, ^ ^, |[ |l^, |ll|; |lj, p^ . 
Another important tool in QCD applications of the chiral ensemble (as well as in studies on 
closely related non-Hermitian random matrices [^) proved to be the following generalisation of 
the Itzykson-Zuber type integral 







1 exp 


4Tr(x. 



VYdU^ 



Here the integration in (||) goes over unitary n x n matrices U,V G U(n) with dii{U, V) being the 
corresponding invariant measures. The matrices Xd and Yd are diagonal: 

Xd = diag{xi,X2,.-.,Xn), Yd = diag (j/i, 2/2, • ■ • , 2/n) (6) 



2 



and /o(a;) denotes the Bessel function of hiiaginary argument and 0^^ order. In fact, the mtegral 
formula (|^) was known aheady to Berezin and Karpelevich , and rediscovered very recently by 
Guhr and Wettig |Q and by Jackson, Sener and Verbaarschot 

A more general correlation function Eq. (Q) can be studied for the particular case of the Gaus- 
sian measure ( chiral Gaussian Unitary Ensemble, eq.(H)) by the procedure known in the literature 
as the supermatrix (or " supersymmetry" ) approach. The method was pioneered by Efetov 
in the theory of disordered systems. First, one represents each of the characteristic polynomials 
in the numerator as the Gaussian integral over anticommuting (Grassmann) variables. Similarly, 
characteristic polynomials in the denominator are represented by the Gaussian integrals over usual 
commuting complex vectors. This allows to average the resulting expressions in a simple way. The 
resulting multiple non-Gaussian integral is brought to a tractable form by the so-called Hubbard- 
Stratonovich transformation. The latter transformation has, in general, a quite complicated ana- 
lytical structure which mixes commuting and anticommuting as well as compact and non-compact 
integration variables. Such a trick allows one to expose the integration variables amenable to the 
saddle-point treatment in the limit of large matrix size N. 

The first applications of the supersymmetry technique to the chiral case ||, showed that 
the computation was quite technically involved even for the lowest correlation functions. In fact, 
any explicit calculation beyond the two-point correlation function of the resolvents seemed to be 
problematic. 



In our paper |g3| we demonstrated that an alternative method suggested in 22 1 enabled us 
to calculate general spectral correlation functions for random Hermitian matrices without chiral 
structure in a systematic and efficient way. Our starting point was the same as for the supersym- 
metry approach. At the same time, we integrated out Grassmann variables at a very early stage 
and avoided the complicated Hubbard-Stratonovich transformation thus seriously departing from 
the general spirit of supersymmetry. 

In the present paper we show that the methods developed in ^ can be successfully applied 
to calculate the general correlation function of products and ratios of characteristic polynomials 
for chiral GUE matrices (for the simplest case rif, = n/ = 1 this fact was already demonstrated in 
11)- 

Our main result is the formula: 

ICn{X) = const ^dct ( / dRe-^"R^-"^+'-^Jkix'-''^^)) (7) 



A{X|}A{X|} 



valid for any values of the parameters N,ni,,nf, Xb = 2NMb, Xp = 2NMf, provided N > nt- 
Here we introduced the matrix X — diag of the size n — Uf + nil and denoted 

' \ Ka{z) ,nf + l<k<n 

where Kq{x) denotes the Macdonald function of 0*^ order. 

Quite remarkably, the final expression is represented in a compact determinantal form reminis- 
cent of that arising from the orthogonal polynomial method for positive moments. In a separate 
publication [p^ we show that such a determinantal form is in no way accidental and is not specific 
for the chiral case. In fact, we managed to derive it in full generality not only for the Gaussian 
distribution, but for arbitrary unitary-invariant potential. For the present case of chiral GUE one 
can perform the asymptotic analysis of Eq.(^ in the "chiral limit" N ^ oo with fixed Xb and Xp 
and obtain more compact determinant expression Eq. (p&) . This is also a new result going beyond 
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the known asymptotic expressions for the product of positive [^0|, |l2| as well as only negative 
moments. 

The key technical tool for our analysis is the matrix integral: 



exp 



dfi{f,f^) = const 



det [Koix^yj]] |i<ij<„ 



A(X?,...,X2) A(2/2,...,y2) 

(8) 
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Now Xd and Yd are two positive definite diagonal matrices: 

X<j = diag(a;i,X2, . . . jXn) > 0, Yd = diag (yi, y2, • • ■ , J/n) > (9) 



The integration in (g) goes over T e GI(n,C)/U(l) x ... x U(l), i.e. over arbitrary n x n com- 
plex matrices with positive diagonal elements, with dfi{T,T'') being the corresponding invariant 
measure. 

The integrals and (||) can be looked at as certain matrix Bessel and Macdonald functions. 
In fact, second integral (||) is a very natural non-compact counterpart of the integral 

The structure of the paper is as follows. First we follow the (improved) version of the method 
suggested in [ p^ and derive a convenient representation of the general correlation function in terms 
of a matrix integral. Then we integrate out irrelevant degrees of freedom exploiting Eqs.(^ and 
(^ and reveal a simple determinantal structure of the resulting expression. As a by-product we 
expose variables amenable to saddle-point treatment in the limit of large matrix size N. Finally 
we calculate the resulting correlation functions in the "chiral scaling" limit N oo, niN — x < oo. 
In the Appendices we present a derivation of the Itzykson-Zuber-like integral (^) using standard 
"diffusion equation" arguments, as well as discuss some technical details of matrix parametrisations 
used in our calculation. 



2 Correlation Functions for Chiral CUE 

As we discussed in the Introduction the starting point of our method is similar to the standard su- 
persymmetry approach [ pT| and is based on a simultaneous exploitation of commuting (or bosonic) 
and anticommuting Grassmann (or fermionic) variables. We introduce 2nf fermionic and 2nb 
bosonic vectors, each of them with N components: 

• Xfe, <^fc - fermionic vectors, fc = 1, 2, . . . , 

• Si, p; - bosonic vectors, / = 1, 2, . . . , n;, 

These vectors enable us to rewrite the ratio of products of characteristic polynomials in Eq. (P as 
the following integral: 

fl2(j,mf) 

!^ = const / VJ'VB exp{~Ajr{J) - Ab{J)} (10) 

riz(j,m«) 

1=1 ^ ' 
where we have denoted: 

"/ "i 

VJ^ =Y\_dxldxkd(pldipk, VB = Yl_dsldsidpldpi (11) 

fe=i 1=1 
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and 



fe=i 



(12) 



(13) 



1=1 



The result of the above manipulation is that we can easily perform the ensemble average over the 
matrices J in Eq.(^) using the identity: 



-Tr{.rC+D.J) 



eirTrCD 



(14) 



where in our case: 

C ^i^si^pl -i^(pk<E)xl > D ^i^pi<»sl -i^Xk'^ifil 

Ik Ik 

After inserting Eq. (|l^) into Eq. (^) and performing the average in the way described above we 
find it convenient to introduce two matrices X,y of the size Uf x Uf with entries: 

Xkik2 = xl,Xk2, 3^fclfe2 = Vli<^fe2> (15) 

as well as two matrices Qbi,Qb2 of the size rib x ni, : 

Q'^'^=slsi,, Q'^^^^pIpi, (16) 
Then the result of the average can be written in a compact form as 

/C(M/, Mb) = const J VTVB exp {-W^ - Wj? - W^g} (17) 

where 

Wjr = TYl^Mf{x + y)-^xy^ (18) 

Wb = Tr I^KUQbi + Qb2) + ^QbiQb2^ (19) 

" ^ 51 ® P?)'^'^- + "^I^P' (20) 

In contrast to the standard supersymmetric approach from now on we deal with fermionic and 
bosonic integration separately. Let us first integrate out fermionic variables in Eq.(17). To re- 
move the non-gaussian terms in the exponent we employ the simplest version of the Hubbard- 
Stratonovich transformation (cf. (|l^): 



exp 



^Trl^i-j)) = const ^ dQi^dQ^, exp {-iVTr(g^QF) - Tr(Qj.i) ~ Tr(Q],3))} (21) 



where Qp \s a. Uf x Uf complex matrix. We insert the above formula to the integral Eq. ( p^ ) and 
change the order of integrations. Then the integration over fermionic vectors can be done explicitly 
yielding the following representation for the correlation function: 



JC{Mf,Mb) = j VBdQpdQ'^p deti exp j-Wg - NT:t[Q'^pQf)^ 



(22) 
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where the {2nfN) x {2njN) matrix A has the foUowing structure: 



A 



(23) 



with rhp = Alf + Qf, and ,B = ^ s; (g) (p[). By simple algebraic manipulations one can show 

1=1 

that: 



det 



Inf Ig 



= dot 



([at«] 



(g) j at (g) Ijv 

for any Uf x Uf matrix d and N x N matrix b. Moreover, it is easy to see that 



N 



det i^Ra ® Iat - 1„ (g i?6 j = n det (r«lAr - = n n 

•i— 1 i—lj — l 



' b 



for any two Hermitian matrices Ra and Rf, with eigenvalues ri^'', ...,ri"^ and r^^-*, ...,r[^\ respec- 
tively. Finally, we notice that the matrix TZb — BB^ has obviously rank rib, i-e. ni, nonzero 
eigenvalues, and verify that TrTZ^ — Tr (^Q biQ B2^ for any positive integer p. As a consequence, 
for Tib < N we have 



det 



In - N-^Ub 1 = 



det (rWl„, -iV-2g^^g^^^ 



and immediately obtain the following formula for the determinant entering the equation Eq.(p2| 



det A = det 







rhpm'p 


det 



det rrVp-mp (g 1„ — 



QbiQ 



B2 



(24) 



We see now that the integrand depends on vectors s/, p; , Z = 1, nb only via the scalar prod- 
ucts S;-^S(2 and PQPi2 forming the matrices (5b(i,2)j see Eq.(^6|). Then one can convert integration 
over vectors into integration over the Hermitian matrices Qbi > 0,Qb2 > using the integration 
theorem, see [E2l ESl: 



VB (Qbi,Qb2) / dQBi / dQB2 \Qbi,Qb2) det (Qbi) det ((3_b2 
J ^ ' JQbi>0 JQb2>0 ^ / I \ / \ 

and to write down the correlation function under consideration as a matrix integral: 



dQ 



Bl 



i>0 



dQF I dQl ei^rx,{QUQ.)-NTT{Q^,Q.) [jet (q^Qf 

dQB2 e-^Trx«(Q-i+Q-)-^Tr(Q«iQ«.) J^^^ (g^^Q^^ 



Qb2>0 



(25) 



X det 



(Q^Qf) g) 1«, - In^ ® (C 



'B1WB2 



where we denoted Xb ~ 2N M\j and Xp — 2N Mf, correspondingly. 
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Next natural step is to use the singular value decomposition Qp = W qV, where U,V ^ U(nf) 
are two different UfXUf unitary matrices, and q = diag(gi, g„j.) are singular values, i.e. positive 
square roots of the eigenvalues of Q'^pQp > 0. The integration measure in the new variables is given 
by dQ^pdQp cx Y\j;Li^k'^^kA^{(f'}d^{tj)d^{V), with dii{U,V) being the corresponding Haar's 
measure on the group U(nf) and A{f} = Uk^<k2 illi -ik)- 

Now we have to find an appropriate parameterisation for the pair of positive definite Hermitian 
matrices Qbi > and Qb2 > of the size ni, x ni,. One can prove that such a pair can always be 
uniquely represented as 



Qbi — 



P 



T 



1B2 



(26) 



in terms of a positive definite diagonal P > and a general complex matrix T with real positive 
diagonal entries: T £ GL(n,C)/U(l) x . . . x U(l). The corresponding integration measure is: 



dQsidQB^ oc \{pidpi J| {pf - pl,fdn{f, ft) , d/x(f , ft) = dfdf^ det 



fft 



(27) 



1=1 



The proof of the above statements is presented in the Appendix A. 

The only terms in the integrand Eq.(p5|) which depend on the matrices U,V and T,T'^ are 
e5TrXF(Qt,+QF) a,nd e^i^^^''('^''^+'^''^\ respectively. It is immediately evident that the in- 
tegration over dfi{U,V) can be performed using the matrix Bessel function Eq.(||) whereas its 
counterpart Eq.(||) allows us to integrate out the matrices f and ft. Moreover, the Vandermonde 
determinant factors arising in the denominator after that integration cancel away one of such fac- 
tors coming from the integration measure. As the result, the integrand becomes proportional to 
fully antisymmetric Vandermonde factor. This antisymmetry can be used to show that the inte- 
gral of the determinant of Bessel functions is equal to the integral of the product of the diagonal 
elements of the corresponding nj x nf matrix, multiplied with the overall n/! factor. (This factor 
counts the number of terms in the determinant). The same procedure is applicable to the integral 
of the rif, X nj, determinant containing Macdonald functions. Denoting = Rb , q^ = Rp with 
eigenvalues r'^^ and i?^-* , respectively, the resulting expression is given by: 



ICn{Xp, Xb) oc 



A{X|}A{X|} 



X / dRpA{Rp} 

'rf>0 



X I dReAlRB} 
'RB>n 



det Rb 



det Rp 



fc=i 



1=1 



e-^TrA. -Q^^ Ui;V4^ n - <) (28) 



k=l 1 = 1 



We further introduce a matrix R — diag(Rp, Rb) of the size n 
X = diag ( Xp, Xb ) of the same size. Then 



rih and another matrix 



)Cn{X) = C 



N 



_ 1 •T'rY2 

e 41V -^'-^F 

A{Xl}A{Xl} Jr>o 



dRA{R} 



detR 



[]jfe(xW/RW) (29) 



fe=i 



where 



Jk[z)-<^ Xo(z) ,nf + l<k<n 
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Finally, the presence of the Vandermonde determinant in Eq. (|29|) suggests to rewrite the corre- 
lation function in the form of n x n determinant. This leads us to the determinant formula Eq.(^ 
which is the principal result of the present paper. Our formula generalises known expressions for 
the product of positive moments | [lO| . It is valid for any values of the parameters N,ni,,nf, Xb, Xp, 
provided N > rib. 

The proportionality constant Cat in the above formulae is fixed by the obvious normalisation 
condition ICNiX)\xB=Xp = 1 following from the very definition (^. In fact, we find it easier 
to restore this constant considering the limit Xb — > 00,^^7 00, when obviously 1Cn{X) —f 

(2jV)2^^(»b-»/) lifc=i . On the other hand, inspecting the integral Eq.(|2|) one notices that 

nr=i i^B \ 

(I) (k) 

in the discussed limit one can effectively put Rg Rp in the integrand. Then the integral 



decouples into the product of two integrals: ICn{X)\xb,Xp^c 
where 



ICBiXB)\xB^ooK.F{XF)\ 



!Cb{Xb) = 



ICf(Xf) 



1 



dRBA{RB} 



B.B>0 



det Rb 



X, 



e 4" 



Trxl 



dRpAlRF} 



B.F>0 



det Rp 



N 



-XTlRf 



fe=i 



(30) 
(31) 



In the limit of large X^' ^ N^^^ one can use the asymptotic expression for the Macdonald function 
Kq{x ^ 1) = in the integrand of (|3^) and safely neglect the terms iVi?^-* in the exponent. 

The resulting integral is of the form 



/•CO i-oo " A/y2 v2l " 

/ .../ A{R,, n R^er-^^dR, = 2" ^ n r [20- 



(32) 



In our case a = N — rth — 1/4 which yields: 



ICb{X 



B) Xb^oo 



(27r)"^/2 H [xjj'^J n r [2(^^ - nfc + - 1/2] 



(33) 



For the second integral we should use the asymptotic expression for the modified Bessel function 
loix > 1) = 



(k) 

. Here, however, we can not neglect the terms NRp in the exponent since it 

(k) 

ensures the convergence of the integral. In fact, it is easy to see that for Xp ^ 2N the integral 
is dominated by small vicinity of the value Rp^ = x'^p^ /2N . Taking into account small Gaussian 
fluctuations around those points we find that 



1Cf{Xf)\xb 



n [^'b\ 



k=l 



2N 



The equations (|3^)-(|3^) yield the following value for the normalisation constant: 



Cn — 



(27r)"^/2 n r [2(iV -Ub + l)- 1/2] 



(34) 
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As usual, we would like to extract the leading behavior of the above correlation function in the 
scaling limit N ^ oo. Here we are mostly interested in the so-called "chiral scaling" /C'^''"'(X) = 
liniTv^oo I^n{X) which amounts to keeping the parameters X fixed when treating the corresponding 
integrals by the saddle-point method. The saddle-points are extremals of the function C{Rk) — 
Rk — lni?fc, k = l,...,n in the domain Rk > 0, which are given by Rk = 1. To pick up a 
nonvanishing contribution one has to take accurately into account Gaussian fluctuations around 
the saddle point. It is done in the most straightforward way by exploiting the expression ( p9[ ) 
and putting there Rk — 1 + £,k, with — oo < < oo parameterizing the fluctuations. Introducing 
S = (Ci, ■■■,£.n) we have: 



1 



A{X|}A{X|,} 
1 



(35) 



A{X2 }A{X|} A^- 



lim det 



fe=i 



Lk=l 



Evaluating the limit we finally find: 



det 



-'0 -^F 



X 



("/) 



Ko Xg 



(1) 



Ko 



X 



("b) 



A{Xl}A{Xl} 



^F ^0 



X 



X 



X 



("/) 



B 







("/ 
F 

(1) 



B 



(^(^".)yn-l^("-l) 



X 



(36) 



where Iq\z) and Kq\z) stand for l-th derivatives of the Bessel and Macdonald functions. 

Such an expression provides us with the most general correlation function for chiral GUE 
and generalises earlier results known for > , rib = |l0|| , Uf = , ni, > |22], and 
nf = ni, = (1, 2)1^. For example, for Uf — ni, — 1 the equation (|36|) amounts to 



IC''''''{Xf,Xb) (X [XFh{XF)Ko{XB) + XbIo{Xf)Ki{Xb)] 



and is simply related to the formula found in Q . 



2.1 Conclusion 

We have demonstrated that the general spectral correlation functions containing both products and 
ratios of the characteristic polynomials of chiral GUE matrices can be calculated in a closed form. 
Our method amounted to representing the characteristic polynomials in terms of the Gaussian 
integrals and exploiting the Itzykson-Zuber type integration formulae Eqs.(|^,^), both over compact 
and non-compact domains. The results were shown to have an attractively simple determinantal 
structure reminiscent of that arising in the method of orthogonal polynomials. In a separate 
publication we demonstrate that indeed there exists a way to calculate the general correlation 
function by a method resorting to the orthogonal polynomial technique. This way makes clear 
that the determinantal structure observed in the present paper has, in fact, a very profound origin 
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and allows one to extend our results to other unitary-invariant ensembles, with or without chiral 
structure. As interesting prospects for future research we would like to mention a challenging 
problems of extending the suggested methods to other symmetry classes (see |5| ) , as well as 
to ensembles of non-Hermitian random matrices. The latter ar important for QCD applications 
(see e.g. Q and references therein) and for the problems of quantum chaotic scattering pTf . 
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3 Appendix A. A parametrisation for the pair of positive 
definite matrices Qbi and Qb2 

Given two positive definite complex Hermitian matrices Qsi > and Qb2 > of the size ni, x ni, 
let us introduce Qbi = [Qbi]~^ > 0, so that 



dQBidQB2 = dQBidQB2 (dct Qbi 



-2n,, 



Now we use that any pair of complex Hermitian Qsi > and Qb2 > can be parameterized 
as (see, for example 

Qbi^TT^ , Qb,^TPT^ 

where the matrix P is a diagonal positive definite: P — diag(pi, ...,pn^) > and the matrix T 
is general complex, with the only restriction that its diagonal entries are real and positive: T S 
GI(nb,C)/U(l) X ... X U(l). The latter condition is necessary to ensure one-to-one correspondence 
between the parametrisations. Our goal is to calculate the related Jacobian. 
For this goal we use the following relation between the matrix differentials: 

dQBi^T[ST + ST^ Tt 

where we denoted ST = T-^dT and ST"^ = dT^ [T^] Analogously 

dQB2 = T [6TP + dP + PST^ 

Further denoting SQbi - P-^dQei [T^^ , SQb2 = T-^dQB2 [T^^' we first calculate the Jaco- 
bian of the transformation from the set of variables [SQbi].^ , [SQB2]ii , [SQBi]i^j , [SQ*Bi]i^j [<5(3b2]j<j , [SQb2\i<j 
to a new set of independent variables dpi, {6T)ii, {6T)i^j, {5T)*^j, {5T)i^j, (5T)*^j where \ < i,j < 
ni, and * stands for the complex conjugation. The calculation is simple and yields the factor that 
we symbolically write as 

S{Qbi,Qb2) YJ/ \2 

We need also the following intermediate Jacobians: 

|g^i#4 = det(rTtf"^ 

0{Qbi,Qb2) 
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det (TT^) 



-(ri6-l/2) 



d{T, Tt) 

Then all the factors taken together yield the full Jacobian: 



diQBi,QB2) _ S{Qbi,Qb2) „ rf(QBi,Qs2) ^ <5(r, Tt) 



n(p.-P,f det(Trt)(-+^/^^ 



<5(gsi,Qs2) d(T,Tt) 



(37) 



i<3 



To arrive to the parametrisation Eq.([2^) used in the main text of the paper we change P —t 
so that dP ~* 2"''dPdetP in the measure and then change T T p-^/'^ so that Tt p-^/'^T^. 
Then Qbi = Qbi = [TP^^T''] = [T^ PP-^ as required and the matrix Qb2 stays equal to 
TPT'^ . The transformation amounts to change in the measure dTdT^ detP-("^-i/2)dTdTt. 
Finally, taking into account the factor 



det 



Q 



Bl 



detp-^"" det (TT^ 



we arrive to the expression for the measure given in Eq.(p7|). Let us note that diJ.{T, T^) is exactly 
the invariant measure on the manifold GI(n,C)/U(l) x ... x U(l). 



Appendix B. Matrix Macdonald functions associated with 
integrals over complex matrices 

The matrix Bessel functions that correspond to Itzykson-Zuber-like integrals over unitary matrices 
were considered in details by Guhr and Wettig, Guhr and Kohler |20j. Below we consider the 
matrix Macdonald functions that are associated with integrals over arbitrary complex matrices 
with real positive diagonal elements. 

Let Xd and Yd be two diagonal matrices with real positive elements, 

Xrf = diag(xi,a;2, . . . ,x„) , diag (j/i, y2, • • • , y™) (38) 

Our aim is to show that the integral: 



HXd,Yd) 







1 exp 


4Tr(x. 



dfi{T) (39) 



(where T are arbitrary n x n complex matrices) can be considered as the matrix Macdonald 
function. We are going to demonstrate the following expression for the function ^{Xd,Yd): 

'^(^^ = A(.?,...,.^)A(,?:..:,,^) 

where Kq{x) is the Bessel function of an imaginary argument of 0^^^ order. 

Given three n x n complex matrices X, A, B we consider the Laplace operator Dx acting on 
complex matrices as: 

Dx^ E (t^h^ + i^^t:?^) (41) 



\d{ReX,jY a(ImX,-^2 
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Then we construct a function Vt^(X, A, B) with the property 

DxW{X, A, B) = Tv{AB)W{X, A, B) 
In particular, the following function 



W{X, A, B) = exp 



]^Tv{XA + BX^) 



(42) 



(43) 



satisfies the Eq.(^) as can be checked by direct calculations. Let us put specifically A = TYdT^ 
and B ^ (Tt)-iyrfr-i in Eq.(|2l). We obtain: 



(44) 



1, 



Dx exp 

Tr(y|)exp 



--Tr {X [TYdT^ + {T^y^YaT-^]) 



-Tr {X [TYdT^ + [T^y^YaT-^]) 



As soon as the integration over complex matrices T commutes with the Laplace operator Dx we 
conclude that the matrix function ^{Xd, Yd) defined by the integral Eq.(|3^) satisfies the following 
differential equation: 

Dx^Xd, Yd) = TY{j2)^{Xd, Yd) (45) 

To derive the explicit formula (Eq.(|39|)) for the matrix function ^(Xd^Yd) we apply the method 
proposed by Guhr and Wettig pfl] . Using the singular value decomposition for an arbitrary complex 
matrix X: 

X = U^XdV, U^U = In, V^V = 1„ (46) 
the radial part Dx^ of the Laplace operator Dx must have the following expression: 

n 

l[xi (47) 



^ ^ i—l 



2 ^2n 

, . . . , 



Guhr and Wettig noted that the radial part Dx^ of the Laplace operator Dx is separable. It 
means that for an arbitrary function f{xi, X2, ■ ■ ■ , Xn) the following identity holds: 



f {x\ , . . . , Xri) 



y{ — 



'' /\{xl,...,xl) A{xl,...,xl) f^^ydxl xkdxkj A{ 



I d \ f{xi, ...,Xn) 



(48) 



The separability of the operator Dx^ enables us to solve the differential equation (Eq.^Sj) and to 
prove formula (^). We use the following ansatz for the solution Yd) of differential equation 

^^ 1'(xi,...,a;n;yi,...,y„) 
<^{Xd, Yd) = ^7-2 o..,,2 



A{xl...,xl)A{yl...,yl) ^^^^ 

We insert the above ansatz to the differential equation (Eq.^Sj). Applying the property (Eq.|4^) of 
the radial part Dx^ of the Laplace operator Dx we find: 

E f^fc + '^^^ "^{xi, ... ,x„;j/i, ...,?/„) = (2/1 +1/2 + •• • +yi) • • • ,a;„;yi, . . . ,y„) (50) 



k=l 
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The form of the above differential equation suggests to fook for a solution as a sum over permuta- 
tions of the set {1,2, ... ,n). Each term in the sum can be taken as a product of n multipliers and 
each multiplier is represented by a (kernel) function taken at different Xi and yj, i.e. 

*(xi, ... ,a;„;yi, ...,?/„) = ^ (j){xi,y^^i)) (j){x2,ya(2)) ■ ■ ■ (f}ixn,yain)) (51) 

o-e s„ 

where Co- is some constant dependent on a particular permutation a. To determine Ca we ex- 
plore the symmetry properties of the function '^{xi, . . . , a:„; yi, . . . , yn). The function ^(Xd, Yd) is 
symmetric under permutations inside the set {yi, y2, . . . , yn} as can be seen from the integral repre- 
sentation (Eq.^). Indeed, the result of any such permutation can be represented by a permutation 
matrix P as 

Yd ^ - PYdP^ (52) 

Let us replace the matrix Yd by Y£ in the integral (Eq.^9|). Changing the variables of integration, 
Ti = TP, we see that $(X<j, PYd-P^) = ^{Xd,Yd), i.e. the integral (Eq.|9|) is symmetric under 
permutations inside the set {yi, y2, . . . , j/n}. As soon as ^{Xd, Yd) and '^{xi, . . . , Xn,yi, ■ ■ ■ , yn) 
are related by Eq.(p9[) the function '^{xi, . . . , Xn', yi, . . . , yn) is antisymmetric under those permu- 
tations. It gives the coefficient Ca in the formula (Eq.|5l|): 

a = i-l)"' (53) 

where ~ 1 for odd permutations and i^o- = for even permutations. We then obtain the following 
expression for the function ^'(xi, . . . , x„; j/i, . . . , y„) from the Eq.(pl|): 

^{Xi, ...,Xn;yi,-- ■,yn) = det {(l){Xi,yj)) \l<i,j<n (54) 

We insert the expression Eq.(|5l]) for '^{xi, . . . , Xn] yi, ■ ■ ■ , yn) to the differential equation (Eq.pO|). 
As a result we find that the kernel (f>{x, y) satisfies the differential equation: 



(t>{x,y) ^y^(t>{x,y) (55) 



Both Bessel functions of an imaginary argument Io{xy) and the Macdonald functions K(){xy) satisfy 
the above differential equation. However the integral representation (Eq.^) dictates the choice 



for the kernel function 4>{x, y). Indeed, we note that for n = 1 the integral (Eq.39) is proportional 
to the kernel function (p{x,y). This integral representation coincides with that for the Macdonald 
function Ko{xy) and we conclude that (f>{x, y) = const Ko{xy). Then equations (|49| ) and (^J) yield 



the desired formula (Eq,40) 
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